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boson or a tt pair. Depending on the production mode, Higgs boson decays to W+W~ , ZZ, bb, 
T^T~, and 77 are examined. We search for a Higgs boson with masses {mH) in the range 90- 
200 GeV/c^. In the absence of a signal, we expect based on combined search sensitivity to exclude 
at the 95% credibility level the mass regions 90 < mH < 94 GeV/c^ 96 < < 106 GeV/c^ 
and 153 < mn < 175 GeV/c?. The observed exclusion regions are 90 < ma < 102 GeV/c^ and 
149 < ruH < 172 GeV/c^. A moderate excess of signal-like events relative to the background 
expectation at the level of 2.0 standard deviations is present in the data for the run = 125 GeV/c^ 
search hypothesis. We also present interpretations of the data within the context of a fermiophobic 
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Within the standard model (SM) [l| of particle physics, 
the mechanism of electroweak symmetry breaking [2] im- 
phes the existence of a single observable particle referred 
to as the Higgs boson, H . The mass of this neutral scalar 
is not predicted by the theoretical framework of the SM 
and must be measured experimentally. Similarly, Yukawa 
couplings between fundamental fermions and the Higgs 
field, which are responsible for fermion masses, are not 
predicted by the SM. 

Precision electroweak measurements from LEP, SLC, 
and the Tevatron have been interpreted within the con- 
text of the SM to constrain the mass of the potential SM 
Higgs boson . Including the most recent W boson and 
top-quark mass measurements from the Tevatron 0, , 
the electroweak data are consistent with a Higgs boson 
mass smaller than 152 GeV/c^ at the 95% confidence 
level, within the framework of the SM. Direct searches 
at LEP exclude the SM Higgs boson for masses less than 
114.4 GeV/c2 [i]. 

Recently, a new particle was observed in data col- 
lected from y^s =7-8 TeV pp collisions at the CERN 
Large Hadron CoUider (LHC) by the ATLAS Q and 
CMS [1| collaborations. The reported measurements of 
the observed particle are consistent with the expecta- 
tions for the SM Higgs boson with a mass of roughly 
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125 GeV/c^. The specific final states contributing the 
greatest amount of significance to these observations are 
77 and ZZ ^ Ml [^]. Complementary evidence was 
also reported recently in the hh final state based on a 
combination of searches from the CDF and DO experi- 
ments [loj . Precision measurements of the properties of 
the new particle such as its spin, parity, production rates 
via the different mechanisms, and decay branching ra- 
tios are necessary for identifying if the new particle is in 
fact the SM Higgs boson. Higgs boson searches at the 
Tevatron obtain most of their sensitivity from produc- 
tion and decay modes that are different from those of the 
LHC searches. Tevatron measurements therefore provide 
important contributions to the available constraints on 
several of these properties. 

The SM Higgs boson production process with the 
largest cross section at the Tevatron is gluon fusion. As- 
sociated production with a. W ot Z boson {VH) is the 
second largest. The cross section for WH production is 
twice that of ZH and is about a factor of ten smaller 
than gluon fusion. The Higgs boson decay branching 
fraction is dominated by iJ — > 66 for the low-mass Higgs 
boson {niH < 135 GeV/c^) and hy H ^ W+W~ for the 
high-mass Higgs boson {iriH > 135 GcV/c^). An inclu- 
sive search for the low-mass Higgs boson in the H ^ hh 
decay channel is extremely challenging because the 66 
production rate through SM processes is many orders of 
magnitude larger than that expected from the Higgs bo- 
son production rate. Requiring the leptonic decay of the 
associated W oi Z boson greatly improves the expected 
signal-to-background ratio in these channels. As a result, 
Higgs associated production followed by the H —>■ bb de- 
cay is the most promising channel in searches for the low- 
mass Higgs boson. For higher-mass Higgs boson searches, 
the H — >■ W^W~ decay, where leptons originate from 
the W boson decays, are the most sensitive. While the 
H bb and H — J- W^W~ search channels provide the 
best sensitivity, searches made in all final states are com- 
bined to obtain the highest possible sensitivity to the SM 
Higgs boson. 

This article presents a combination of CDF searches 
for the SM Higgs boson. The combined searches incor- 
porate potential contributions from Higgs boson produc- 
tion via gluon fusion, production in association with a 
W ox Z boson, vector-boson fusion (VBF) production, 
and production in association with a top-quark pair. 
Higgs boson decay modes considered are iJ — > W^W~ ^ 
H ZZ, H ^ hb, H ^ T+T-, and H 77. The 
individual searches are performed for potential Higgs bo- 
son masses in the range from 90 to 200 GcV/c^ us- 
ing non-overlapping data sets defined by distinct final 
states. For each search sub-channel, SM backgrounds 
are estimated and validated using data events populat- 
ing appropriately-defined control regions. Finally, a dis- 
criminant, which is typically the output of a multivariate 
algorithm constructed from kinematic event variables, is 
used to separate a potential signal from much larger back- 
ground event contributions. The multivariate algorithms 



are separately optimized for each Higgs boson mass hy- 
pothesis and for each analysis sub-channel. Search re- 
sults are combined by constructing a combined likelihood 
function based on the final discriminant distributions in 
each search sub-channel, taking into account the corre- 
lations among channels. After performing the combined 
search over the full Higgs boson mass range, we focus on 
the 125 GcV/c^ mass hypothesis, motivated by the re- 
cent ATLAS and CMS observations 0, [1] . Assuming the 
LHC signal is present in CDF data, we constrain fermion 
and boson couplings to this new particle. 

We additionally interpret the search results within the 
context of a fermiophobic Higgs model (FHM) (TT| - [l3 |. 
which assumes SM couplings to the Higgs boson except 
in the case of fermions, for which the couplings are as- 
sumed to vanish. In this model, gluon-fusion produc- 
tion is highly suppressed, while branching fractions for 
H — >• 77, H — >■ W^W~ , and H — >• ZZ are enhanced. 
We also consider an extension of the SM incorporating 
a heavy fourth generation of fermions (SM4). Within 
this model, gluon-fusion production is significantly en- 
hanced 

This article is organized as follows: Section briefly 
describes the CDF II detector and the data samples used 
for this combination; Section IIIII describes the predic- 
tions for Higgs boson production and decay that are as- 
sumed throughout, as well as the Monte Carlo models 
used to predict the differential distributions; Section ITVl 
describes the search channels included in the combina- 
tion; Section |V] describes the dominant sources of un- 
certainty in each channel and the correlations of these 
uncertainties between channels; Section IVll describes the 
statistical methods used; Section rVIII presents results in 
the context of the SM; Section lVlIII prescnts results in the 
context of the fermiophobic model; Section IIXI presents 
results in the context of the SM4 model; Section |X] de- 
scribes the measurement of fermion and boson couplings 
in the context of a new 125 GeV/c^ boson; and Section lXj 
summarizes the article. 



II. THE CDF II DETECTOR AND THE FULL 
CDF DATA SET 

The CDF II detector is described in detail else- 
where [H, [1^ . Silicon-strip tracking detectors [l^] sur- 
round the interaction region and provide precision mea- 
surements of charged-particle trajectories in the range 
1 77 1 < 2 [2l|. A cylindrical drift chamber provides full 
coverage over the range I77I < 1. The tracking detectors 
are located within a 1.4 T superconducting solcnoidal 
magnet with field oriented along the beam direction. 
The energies of individual particles and particle jets are 
measured in segmented electromagnetic and hadronic 
calorimeter modules arranged in a projective tower geom- 
etry surrounding the solenoid. Ionization chambers are 
located outside of the calorimeters to help identify muon 
candidates [22| . The Tevatron collider luminosity is mca- 
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surcd with multi-cell gas Cherenkov detectors [23|. The 
total uncertainty on luminosity measurements is ±6.0%, 
of which 4.4% originates from detector acceptance un- 
certainties and 4.0% is due to the uncertainty on the 
inelastic pp cross section . 

All of the results combined in this article, with the ex- 
ception of the H — )- r+T" search, use the full CDF Run II 
data sample. Small variations in the luminosities re- 
ported for the different search channels reflect the appli- 
cation of channel-specific data-quality criteria designed 
to ensure proper data modeling. For example, the silicon 
detector is required to be operational for the H ^ bb 
searches, for which the identification of secondary track 
vertices from b hadrons plays an important role, but not 
in the case of the H ^ 'yy search. The 7? — >■ 77 search 
makes use of the largest data set, 10 fb"^, which is about 
82 % of the 12 fb^^ that was delivered by the Tevatron 
collider and about 99.5 % of the luminosity in which the 
CDF detector was considered to be operational. 

CDF uses a three-level online event selection system 
(trigger) to select beam collision events at a rate that 
can be written into permanent storage. The first trigger 
level relies on special-purpose hardware [2^ to reduce 
the effective beam-crossing frequency of 1.7 MHz to an 
event rate of approximately 15 kHz. The second level 
uses a mixture of dedicated hardware and fast software 
algorithms to further reduce the event rate to roughly 
1 kHz. Events satisfying level-two trigger requirements 
are read out of the detector and passed to an array of 
computers running fast versions of offline reconstruction 
algorithms, which allow for third-level trigger decisions 
based on quantities that are nearly the same as those 
used in offline analyses [2^. The final rate of events 
written into permanent storage is approximately 100 Hz. 
The basic trigger criteria for events used in these searches 
are the presence of high-transverse momentum (pt) [HI 
leptons, clustered calorimeter energy deposits associated 
with partons originating from the collision (jets) |27|. and 
large imbalances in the transverse energies (Et) [21[ of 
measured depositions within the calorimeter, associated 
with evidence for undetected neutrinos within the event 



III. STANDARD MODEL HIGGS BOSON 
SIGNAL PREDICTIONS 

In order to conduct the most sensitive Higgs boson 
search possible, we include contributions from all signifi- 
cant production modes that are expected to occur at the 
Tevatron. When conducting the search using multiple 
production modes, the predictions of the relative con- 
tributions of each mode and the uncertainties on those 
predictions are required. In addition, because we use 
multivariate analysis techniques, the predictions of the 
kinematic distributions for the signal model are also im- 
portant. Here we provide a summary of the tools we 
use to make predictions for the Higgs boson signal. The 



theoretical uncertainties on the signal model play a sig- 
nificant role at higher masses where gluon fusion is the 
major production mode, but are less important for low- 
mass Higgs boson searches where associated production 
is most important. 

To predict the kinematic distributions of Higgs boson 
signal events, we use the pythia Monte Carlo pro- 
gram, with CTEQ5L [2^ leading-order (LO) parton dis- 
tribution functions (PDFs) . We scale these Monte Carlo 
predictions to the highest-order cross section calculations 
available. The pythia differential distributions for some 
important variables, such as the Higgs boson px and 
the number of associated jets, are also corrected based 
on higher-order calculations as described below. The 
gg — > H production cross section is calculated at next- 
to-next-to leading order (NNLO) in quantum chromody- 
namics (QCD) with a next-to-next-to leading log (NNLL) 
resummation of soft gluons; the calculation also includes 
two-loop electroweak effects and handling of the running 
6-quark mas s fsol l3l| . The numerical values in Table |T] 
are updates [33~of the predictions in [3l] | with nit set to 
173.1 GeV/c^ 0, and with a treatment of the massive 
top and bottom loop corrections up to next-to-leading- 
order (NLO) and next-to-leading-log (NLL) accuracy. 
For these calculations the factorization scale (/if) and 
renormalization scale (p-r) are set to fip = fiR = mn, 
and the MSTW 2008 NNLO PDF set [IIL, as recom- 
mended by the PDF4LHC working group [SJ], is used. 
The calculations are refinements of earlier NNLO calcu- 
lations of the gg ^ H production cross section 35 
Electroweak corrections were computed in Refs. [38, 



Soft gluon resummation was introduced in the prediction 
of the gg H production cross section in Ref. [i^ . The 
gg H production cross section depends strongly on the 
gluon PDF and the value of the strong interaction cou- 
pling constant corresponding to the value q of transferred 
momentum, aglji^)- 

For analyses that consider inclusive gg — H produc- 
tion but do not treat events separately based on the num- 
ber of reconstructed jets, we use inclusive calculation un- 
certainties from simultaneous variation of the factoriza- 
tion and renormalization scales by factors of two. We use 
the prescription of the PDF4LHC working group [s^] for 
evaluating PDF uncertainties on the inclusive produc- 
tion cross section. QCD scale uncertainties that affect 
the cross section through their impacts on the PDFs are 
included as a correlated part of the total scale uncer- 
tainty. The remainder of the PDF uncertainty is treated 
as uncorrelated with the QCD scale uncertainty. 

For analyses seeking gg ^ H production that di- 
vide events into categories based on the number of re- 
constructed jets, we follow Refs. [4l|, |42j for evaluating 
the impacts of the scale uncertainties. We treat the 
QCD scale uncertainties obtained from the NNLL inclu- 
sive [30l . l3]|. N LO one or more jets and NLO two or 
more jets [4J| cross section calculations as uncorrelated 
with one another. We then obtain QCD scale uncertain- 
ties for the exclusive gg H + jet, 1 jet, and 2 or more 
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jet categories by propagating the uncertainties on the 
inclusive cross section predictions through the subtrac- 
tions needed to predict the exclusive rates. For example, 
the H+0 jet cross section is obtained by subtracting the 
cross section for production of Higgs bosons with one or 
more jets at NLO from the inclusive NNLL+NNLO cross 
section. We assign three separate, uncorrclatcd scale un- 
certainties with correlated and anticorrclated contribu- 
tions among exclusive jet categories. The procedure in 
Ref. is used to determine PDF model uncertainties. 
These are obtained separately for each bin of jet multi- 
plicity and treated as 100% correlated among jet bins. 

The scale choice affects the pr spectrum of the Higgs 
boson when produced in gluon-gluon fusion, thus bias- 
ing the acceptance of the selection requirements and also 
the shapes of the distributions of the final discriminants. 
The effect of the acceptance change is included in the cal- 
culations of Ref. [1^ and Ref. [4J] , as the experimental 
requirements are simulated in these calculations. The ef- 
fects on the final discriminant shapes are obtained by 
reweighting the pt spectrum of the Higgs boson pro- 
duction in the Monte Carlo simulations to higher-order 
calculations. The Monte Carlo signal simulation used 
is provided by the LO generator pythia which in- 
cludes a parton shower and fragmentation and hadroniza- 
tion models. We reweight the Higgs boson pT spectra in 
the PYTHIA Monte Carlo samples to that predicted by 
HQT [45| when making predictions of differential distri- 
butions oi gg ^ H signal events. To evaluate the impact 
of the scale uncertainty on the differential spectra, we use 
the RESBOS [H] generator, apply the scale-dependent dif- 
ferences in the Higgs boson pT spectrum to the HQT pre- 
diction, and propagate these to the final discriminants as 
a systematic uncertainty on the shape, which is included 
in the calculation of the limits. 

We include all significant Higgs boson production 
modes in the H W+W^, H ZZ, and iJ 77 
searches. Besides gluon-gluon fusion through virtual 
quark loops [ggH), we include Higgs boson production 
in association with slWotZ vector boson (VH) or with 
a top-antitop quark pair (ttH), and vector boson fusion 
(VBF). For the H ^ bb searches, we target the WH, ZH, 
VBF, and UH production modes with specific searches. 
In addition to the leading signal component in each fi- 
nal state, we include contributions of sub-leading signal 
components, which lead to increased signal acceptance. 
The predictions for the WH and ZH cross sections are 
taken from Ref. [13]. This calculation starts with the 
NLO calculation of v2hv jH and includes NNLO QCD 
contributions [4^ , as well as one-loop elcctroweak correc- 
tions [53|. A similar calculation of the cross section is 
available in Ref. [5lj . The VBF cross section is computed 
at NNLO in QCD in Ref. [52|. Electroweak corrections 
to the VBF production cross section are computed with 
the HAWK program |53|, and are small and negative (2- 
3%) in the Higgs boson mass range considered here. We 
include these corrections in the VBF cross sections used 
for this result. The ttH production cross sections we use 



are from Ref. |54{ . 

We use the predictions for the branching ratios of the 
Higgs boson decay from Refs. [H, [s^]. In this calcula- 
tion, the partial decay widths for all Higgs boson decays 
except to pairs of W and Z bosons arc computed with 
HDECAY [5^, and the W and Z pair decay widths are 
computed with Prophecy4f [s^]. The relevant decay 
branching ratios are listed in Table HI The uncertainties 
on the predicted branching ratios from uncertainties in 
the charm- and bottom-quark masses, as. and missing 
higher-order effects are presented in Refs. fssl. [soj. 



IV. SEARCH CHANNELS 

Individual searches typically consist of an event selec- 
tion based on the topology and kinematic properties of 
the final state for the specific Higgs boson production and 
decay mode under consideration. Separation of a poten- 
tial signal from the remaining background contributions 
is obtained in most cases by performing a fit, using the 
signal and background models, for a single discriminant 
variable that is the output of a multivariate algorithm, 
which considers many kinematic event variables as its in- 
puts. The quality of the background model prediction 
for the distribution of each input variable and the fi- 
nal discriminant is studied using orthogonal data samples 
carefully selected to validate the modeling of the major 
background components of each analysis channel. The 
search samples of each analysis are divided into various 
sub-channels based on event information such as types 
of reconstructed leptons, jet multiplicity, and 6-quark- 
tagging characterization criteria. A summary of the in- 
dividual searches and the sub-channels included within 
each is given in Table [Hi We attempt to group events 
with similar signal-to-background ratios within individ- 
ual sub-channels to optimize search sensitivities. This ap- 
proach allows the inclusion of information from less sen- 
sitive event topologies without degrading the overall sen- 
sitivity (for example, events containing higher-impurity 
lepton types). In addition, the isolation of specific sig- 
nal and background components within individual sub- 
channels allows for further optimization of the multivari- 
ate discriminants trained for each, leading to additional 
gains in search sensitivity. The final multivariate dis- 
criminants are separately optimized for each Higgs boson 
mass hypothesis in 5 GeV/c^ steps within the mass range 
under consideration. 



A. H ^ hb searches 

For searches focusing on the H ^ bb decay, the effi- 
ciency for identifying reconstructed jets originating from 
b quarks and the resolution of the invariant-mass mea- 
surement from the two 6-quark jets are of primary im- 
portance. The three most sensitive searches in this decay 
mode utilize a recently-developed multivariate 6-quark- 



TABLE I: The production cross sections and decay branching fractions for the SM Higgs boson assumed for the combination. 
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TABLE II: Luminosities, explored mass ranges, and references for the various processes and final states {i represents e or ^ 
and Thad denotes a hadronic tau-lepton decay) for combined analyses. The generic labels "Ix", "2x", "3x", and "4x" refer to 
separations based on lepton or photon categories. 



Channel 


Luminosity 


juh range 


Reference 




(fb-i) 


(GeV/c^) 




WH — > lubb 2-jet channels 4x(5 b-tag categories) 


9.45 


90-150 


[60] 


WH — >■ lubb 3-jet channels 3x(2 &-tag categories) 


9.45 


90-150 


[60] 


ZH — >■ vvbb (3 6-tag categories) 


9.45 


90-150 


m 


ZH — ;> t^l'bb 2-jet channels 2x(4 6-tag categories) 


9.45 


90-150 


[62] 


ZH — > £'^£~bb 3-jet channels 2x(4 &-tag categories) 


9.45 


90-150 


[62] 


H -> W+W 2x(0 jets)+2x(l jet) + lx(2 or more jets)-f Ix (low-m«) 


9.7 


110-200 


[63] 


H ^ W+W- (e-rh.d) + (/^-rhad) 


9.7 


130-200 


[63] 


WH — >• WW^W' (same-sign leptons)-f (tri-leptons) 


9.7 


110-200 


[63] 


WH WW+W- (tri-leptons with 1 Thad) 


9.7 


110-200 


[63] 


ZH — >■ ZW^W~ (tri-leptons with 1 jet)-|-(tri-leptons with 2 or more jets) 


9.7 


110-200 


[63] 


H ^ ZZ (4 leptons) 


9.7 


120-200 


M 


H — s> T^T^ (1 jet) + (2 or more jets) 


6.0 


100-150 


[65] 


WH + ZH jjbb (2 b-tag categories) 


9.45 


100-150 


[66] 


H ^ '•/'•/ lx(0 jet) + lx(l or more jets) -1-3 x (all jets) 


10.0 


100-150 


[67] 


ttH — >• WWbbbb (4 jet, 5 jet, > 6 jet)x(5 6-tag categories) 


9.45 


100-150 


[68] 
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tagging algorithm (HOBIT) [6^ which is based on the 
kinematic properties of tracks associated with displaced 
decay vertices and other characteristics of reconstructed 
jets sensitive to the flavor of the initiating parton. Tight 
and loose operating points are defined for this algorithm. 
For example, the loose operating point is found to have a 
6-quark tagging efficiency of approximately 70% and an 
associated misidentification rate for light quarks and glu- 
ons of approximately 5%. Compared to the SECVTX 
algorithm, the most commonly used 5-quark tagging al- 
gorithm at CDF, the new algorithm improves 6-tag effi- 
ciency by roughly 20%, for operating points with equiv- 
alent misidentification (mistag) rates. The secondary 
channels that require 6-jet identification (the all-hadronic 
and tin — > ttbb searches) were not updated to use the 
HOBIT tagger and instead rely on the the SECVTX and 
JETPROB [zOl algorithms. The decay width of the Higgs 
boson is expected to be much smaller than the experi- 
mental dijet mass resolution, which is typically 15% of 
the mean reconstructed mass. The H ^ bb searches 
are most sensitive in final states that include two jets. 
However, sometimes initial-state or final-state radiation 
can produce a third jet in the event. Including three-jet 
events increases signal acceptance and adds sensitivity, 
motivating the inclusion of these events in the H ^ bb 
searches. Since a SM Higgs boson signal would appear 
as a broad enhancement in the reconstructed mass dis- 
tribution of candidate 6-quark-jet pairs, dedicated efforts 
to improve the mass resolution have been performed in 
each sub-channel (7l| . Along with improved fe-jct identi- 
fication and jet-energy resolution, the primary iJ — > 66 
analyses have all benefited from increased trigger accep- 
tance by including events from many different trigger 
paths. In many cases, the complicated correlations be- 
tween kinematic variables used in the trigger decision are 
modeled with a neural network using linear regression 
based on event kinematic properties and geometric ac- 
ceptance [72j . 



1. WH lubb search 

The search focusing on the WH — >• ivbb production 
and decay mode [g^l has separate analysis channels for 
events with two and three reconstructed jets. Events are 
further separated into sub-channels based on the type 
of reconstructed lepton and the quality of the tagging 
information associated with the candidate 6-quark jets. 
In particular, separate sub-channels are used for events 
containing a high-quality central muon or central electron 
candidate, a forward muon candidate, a forward electron 
candidate, and a looser central electron or muon candi- 
date based on the presence of an isolated track [t^ [tI] . 
The final two-lepton categories, which provide some ac- 
ceptance for lower-quality electrons and single prong tau 
decays, are considered only in the case of two-jet events. 
For two-jet events, five sub-channels arc used associated 
with each lepton category based on the quality of the 6- 



tagging information associated with each jet: two tight 
tags (TT), one tight and one loose tag (TL), a single 
tight tag (Tx), two loose tags (LL), and a single loose tag 
(Lx). In the case of three-jet events, only two 6-tag sub- 
channels, TT and TL, are considered since the other cate- 
gories contribute negligibly to the overall search sensitiv- 
ity. A Bayesian neural network is used to distinguish po- 
tential Higgs boson signal events from other background 
contributions. 



2. ZH £+rbb search 

The search for ZH — > £^£^bb production and de- 
cay [g^ is based on events with two isolated leptons and a 
minimum of two jets. A combination of triggers based on 
electromagnetic energy clusters and signals in the muon 
chambers matched to reconstructed tracks are used to 
to select events containing Z — >■ ee and Z — >■ /x^ candi- 
dates. Some triggers based on missing transverse energy 
requirements arc also used to select Z — > fifi candidates, 
taking advantage of the apparent imbalance in transverse 
energies that results from the muons depositing only a 
small fraction of their energies in the calorimeter. Neu- 
ral networks are used to select di-electron and di-muon 
candidates [zH . The absence of missing energy from neu- 
trinos allows for improved dijet mass resolution through 
event- wide transverse momentum constraints. These are 
incorporated through corrections to the measured jet en- 
ergies based on the observed 1^^ using a neural-network 
approach. The search maintains separate analysis chan- 
nels for events with two and three jets, as well as for 
events with di-electron and di-muon candidates. Each of 
the resulting four channels is further split into four sub- 
channels (TT, TL, Tx, and LL) based on the quality of 
the tagging information available from the multivariate 
algorithm. Signal is separated from background in multi- 
ple phases. First, three networks are used to distinguish 
the ZH signal from each of the ttb, Z-|-jets, and diboson 
backgrounds. Then, a final network further separates the 
signal from all backgrounds simultaneously 

3. WH, ZH -> ^rbb search 

The search for WH, ZH — >■ ^^^66 production and de- 
cay [6l|, [zl] is based on events with large J^^ and no 
isolated-lepton candidates. Additional background sup- 
pression techniques arc applied to reduce large back- 
ground contributions from multi-jet production pro- 
cesses. In particular, prior to construction of the final 
discriminant, a requirement on a multivariate discrimi- 
nant trained specifically for separating potential signal 
from the multi-jet background is applied to the event 
sample. Events that do not satisfy this requirement are 
used to normalize the remaining multi-jet background 
contribution, which is modeled using a mistag rate func- 
tion for gluon and light-quark jets measured in data and 
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applied to the untagged jets in data events that other- 
wise satisfy the kinematic selection criteria. One of the 
inputs to this multivariate discriminant is a track-based 
missing-transverse-momentum calculation that discrimi- 
nates against false [t^]. A second, final multivariate 
discriminant is used to separate the potential Higgs bo- 
son signal from the remaining backgrounds, such as 4- 
hcavy flavor jets (where heavy flavor refers to jets orig- 
inating from & or c quarks) and ti production. Events 
with two and three jets are treated as a single search 
channel that is subdivided into three sub-channels (TT, 
TL, and Tx) based on the quality of tagging information 
from the multivariate algorithm. 

4. All-hadronic search 

The all-hadronic search focuses on WH,ZH and 
VBF production contributing to the jjbb final state. We 
use events containing four or five reconstructed jets, at 
least two of which have been tagged as 6-quark can- 
didates based on information from the previously de- 
veloped SECVTX [Hj and JETPROB algorithms. 
The use of these two algorithms results in two search 
sub-channels containing events with cither two SECVTX 
tagged jets (SS) or one SECVTX tagged jet and one JET- 
PROB tagged jet (SJ). Large multi-jet background con- 
tributions are modeled from the data by applying a mea- 
sured mistag probability to the non-6-tagged jets within 
data events that contain a single ^-tagged jet but oth- 
erwise satisfy event selection requirements. The multi- 
variate discriminants used to separate potential signal 
from the large background contributions are based on 
kinematic variable inputs including some variables de- 
veloped to distinguish the reconstructed jets originating 
of b quarks from those of light quarks and gluons. 

5. ttH — > ttbb search 

The search for ttH — >• ttbb production and decay 
is based on events with one reconstructed lepton, large 
and four or more reconstructed jets in which at least 
two jets are identified as 6-quark candidates based on 
the SECVTX [H or JETPROB (zOj algorithms. Events 
containing four, five, and six or more jets are analyzed 
as separate channels, and the events within each channel 
are further separated into five sub-channels (SSS, SSJ, 
SJJ, SS, and SJ), based on the number of tagged jets 
and the algorithms contributing to each tag. Multivari- 
ate discriminant variables are used to separate potential 
signal from the dominant tt background contributions. 

B. H ^ T^T~ search 

The search for Higgs bosons decaying to tau lep- 
ton pairs [65| incorporates potential contributions from 



all four production modes. The search is based on 
events containing one electron or muon candidate and 
one hadronically-decaying tau-lepton candidate. To help 
reduce Z/^* — >■ r+r" background contributions, events 
are also required to contain at least one reconstructed 
jet. Events that contain one jet and two or more jets are 
treated as independent search sub-channels. Boosted de- 
cision trees are trained for both sub-channels to separate 
potential signal events from those associated with each 
significantly contributing background production mech- 
anism. Significant numbers of background events are re- 
moved from the samples by placing lower cuts on the 
outputs of each boosted decision tree. The output of the 
boosted decision tree trained to separate potential signal 
from Z/^* — > r+r~ background contributions is used as 
the final discriminating variable for events surviving all 
of the selection criteria. 



C. H ^ W+W- search 

In the search for Higgs bosons decaying to W boson 
pairs [i^l the greatest sensitivity originates from Higgs 
bosons produced through gluon fusion; however, the sig- 
nal contributions from all four production modes are in- 
cluded. The primary search is based on events with two 
oppositely-charged isolated leptons and large -^y, focus- 
ing on the H W^W" — )- i'^vi^v decay mode. The 
presence of neutrinos in the final state prevents an accu- 
rate reconstruction of the candidate Higgs boson mass, 
and separation of a potential signal from background con- 
tributions is based on other kinematic variables. In par- 
ticular, the distribution of angular separations between 
the final-state leptons produced in the decays of W~^W~ 
pairs is significantly different for pairs originating from 
a spin-zero particle, such as the Higgs boson signal, and 
the major backgrounds. 

Events in the primary search are separated into eight 
sub-channels based on the types of reconstructed leptons, 
the number of reconstructed jets, and the invariant mass 
of the dilepton pair. In the case of events with two elec- 
tron or muon candidates, separate analysis channels are 
used for those events with zero, one, and two or more 
reconstructed jets. This separation helps to isolate po- 
tential signal contributions associated with the four sig- 
nal production mechanisms as well as specific background 
contributions such as ti production, which is dominant 
for events containing two or more jets. Based on this 
separation, the final multivariate discriminant used for 
each channel is optimized, leading to a significant im- 
provement in the overall search sensitivity. 

In the case of events with zero or one reconstructed jet, 
separate search sub-channels are used for events contain- 
ing low-purity and high-purity lepton types. Events con- 
taining forward electron candidates, for example, have 
much higher background contributions from W^-|-jets and 
W-\-j production processes where a jet or photon mimics 
the signature of an isolatcd-lepton candidate. 
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A separate search sub-channel is used for events in 
which the dilepton mass is smaUer than 16 GeV/c^. 
The main background event contribution in this kine- 
matic region originates from production, and ad- 
ditional search sensitivity is obtained from the use of a 
separately-trained multivariate discriminant focused on 
separating the potential signal from this particular back- 
ground. Two additional search sub-channels are used 
for events with one electron or muon candidate and 
a second, oppositely-charged hadronically-decaying tau- 
lepton candidate. These event samples contain signifi- 
cant background contributions from VF-|-jets and multi- 
jet production processes, necessitating the use of inde- 
pendent search channels. No further separation of events 
based on the number of reconstructed jets is performed 
within these additional sub-channels. 

Higgs boson production in association with a. W ot Z 
boson in conjunction with the decay H — )- W~^W~ leads 
to additional potential signal events in other, more ex- 
otic final states. The signal contributions are expected 
to be small, but these final states contain much smaller 
contributions from SM background processes. Hence, 
the inclusion of these additional sub-channels improves 
the overall search sensitivity. A search for W'^H — >■ 
W^W~^W~ — ?> i^i/i^i/jj production and decay is in- 
cluded through a sub-channel focused on events contain- 
ing two same-sign, isolated lepton candidates and one 
or more reconstructed jets. Two additional sub-channels 
are used to search for even smaller potential signal contri- 
butions from W+H W+W+W -> l+vi+vi'v pro- 
duction and decay. These sub-channels, one based en- 
tirely on electron and muon candidates and the other 
requiring exactly one hadronically-decaying tau-lepton 
candidate, focus on events with a total of three isolated 
lepton candidates. In all three sub-channels, the final 
multivariate discriminants for separating potential signal 
from other background contributions incorporate multi- 
ple kinematic event variables including the observed 
The provides provides good separation against domi- 
nant background contributions with misidentified lepton 
candidates because of the presence of multiple neutrinos 
within each signal final state. 

Finally, we use events with three isolated lepton candi- 
dates to search for ZH ZW+W l+(.-l+vjj pro- 
duction and decay. Three-lepton events that are found 
to contain a same-flavor, opposite-sign lepton pair with 
a reconstructed mass within 10 GeV/c^ of the Z boson 
mass are classified into one of two separate sub-channels 
based on the presence of one reconstructed jet or two or 
more reconstructed jets. Within the second sub-channel, 
all final-state particles from the Higgs boson production 
and decay are reconstructed (the transverse momentum 
components of the neutrino are obtained from the ob- 
served $rp) and a reconstructed Higgs boson mass is used 
as one of the kinematic input variables to the final mul- 
tivariate discriminant. 



D. H 



ZZ search 



The search for Higgs bosons decaying to Z boson 
pairs l6J] is based on events with four reconstructed 
lepton candidates (electrons or muons). The selected 
events consist primarily of the background from non- 
resonant diboson production of Z*/Z-boson pairs. A 
four-lepton invariant mass discriminant is used for sep- 
arating the potential Higgs boson signal from the non- 
resonant ZZ background. The event Ipj, is used as an 
additional discriminating variable to improve sensitivity 
to potential four-lepton event signal contributions from 
ZH ZW+W- £+e-£+iy£-i^ and ZH e+i-T+r' 
production and decay. 



E. 



— s> 77 search 



The search for Higgs bosons decaying to photon 
pairs (67| incorporates potential signal contributions from 
all four Higgs production mechanisms. Photon can- 
didates are reconstructed in both the central and for- 
ward calorimeters. Conversion (7 — s- e'^e~) candidates 
are also reconstructed in the central calorimeter. Four 
search channels based on these candidate types (central- 
central, central-forward, central-conversion, and forward- 
conversion) arc formed from the inclusive diphoton event 
sample. In order to better optimize the most sensitive 
search category, central-central events are further sepa- 
rated into two sub-channels consisting of events with zero 
reconstructed jets (where the majority of ggH events are 
expected) and one or more reconstructed jets (where the 
majority of VH and VBF events are expected). For these 
sub-channels, multivariate discriminants using the recon- 
structed diphoton mass and other kinematic event vari- 
ables as inputs arc used to separate the potential signal 
from the non-rcsonant backgrounds. In the other analy- 
sis channels, the diphoton invariant mass is used as the 
sole kinematic discriminant. For each Higgs boson mass 
hypothesis, the signal region is defined to be at least ±2 
standard deviations of the expected Higgs boson dipho- 
ton mass resolution. The width of signal windows were 
taken to be 12 GeV/c^ 16 GeV/c^, and 20 GeV/c^ for 
mass hypotheses of 100-115 GeV/c^, 120-135 GeV/c^, 
and 140-150 GeV/c^, respectively. The sideband regions 
around each signal search window are used to normal- 
ize background contributions within the signal region for 
all sub-channels, and to validate the background model- 
ing of the multivariate discriminants for central-central 
events. 



V. SYSTEMATIC UNCERTAINTIES 

The Higgs boson signal production rate is expected to 
be small compared with the copious backgrounds pro- 
duced in pp collisions at Tevatron energies. Systematic 
uncertainties associated with background predictions can 
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be significant relative to expected signal rates in the high- 
est s/b bins of the discriminant distributions. Therefore, 
it is expected that systematic uncertainties can have a 
large impact on search sensitivity. As an example, in the 
case of the search for a 125 GeV/c^ Higgs boson the inclu- 
sion of systematic uncertainties weakens the sensitivity of 
the combined analysis by roughly 20%. We consider un- 
certainties that affect the normalizations as well as those 
that affect the shapes of the multivariate discriminants 
used in the searches. We refer to these respectively as rate 
and shape uncertainties. Some systematic uncertainties 
are correlated between analyses, between sub-channels 
within an analysis, and between signal and background 
predictions within a sub-channel. The nature of the fits 
that are performed requires careful evaluation of the com- 
mon and independent sources of systematic uncertainty. 
The details of the statistical treatment of the uncertain- 
ties are described in Sec. IVII 

The most important rate uncertainties in the back- 
grounds to the WH Whh and ZH Zbb searches 
come from the VF-|-jets and Z-|-jets backgrounds. These 
uncertainties are separated into heavy-flavor components 
and mistags. The mistags are calibrated using data con- 
trol samples. Because these backgrounds are calibrated 
in situ in events with different selection requirements 
from the analysis search region, and because the differ- 
ences between the predictions and the true rates may not 
be the same between the W^-|-jets and Z-|-jets samples, 
the M^-l-heavy flavor and Z+hea,vy flavor uncertainties 
are not correlated between analyses, but are correlated 
between sub-channels of a single analysis. This treat- 
ment ensures that we do not use the .Z-|-jets searches to 
cross-calibrate the backgrounds in the VK-|-jets searches 
and vice versa. 

The uncertainties on the 6-tag efficiencies for each 6-jet 
selection requirement are evaluated both for true b jets 
and for mistagged jets. These uncertainties are propa- 
gated to each 6-tag category. The resulting uncertainties 
are treated as correlated between the signal predictions 
and the background predictions. The uncertainties re- 
lated to the 6-tag efficiencies are treated as correlated 
between analyses that use the same 6-tag algorithm. Sim- 
ilarly, the uncertainty on the total integrated luminosity 
as measured by the luminosity monitor is considered cor- 
related among all signal and background MC-bascd pre- 
dictions in all analyses. 

We ensure that each analysis uses the same cross sec- 
tion assumptions and theoretical uncertainties on the 
prediction for the production of diboson [t^] 
W^Z, and ZZ), tt [tI], and s-channel and t-channel 
single-top quark [t^ events. The three diboson pro- 
cesses share common dependencies on factorization and 
renormalization scales and PDFs, so we correlate the un- 
certainties on all three production modes, and correlate 
these uncertainties across all channels. 

The jet energy scale is calibrated with experimental 
data using events in which a photon recoils from a jet, 
and events in which a leptonically-decaying Z boson re- 



coils from a jet (27|. The associated uncertainties are ap- 
plied to each analysis. They change the predicted rates 
of events passing the respective selections, largely due 
to jet Et requirements, but also distort the predicted 
shapes of the distributions of the final discriminant vari- 
ables. Hence, the systematic uncertainty from the jet 
energy scale can be further constrained in situ. We do 
not, however, correlate the jet energy scale uncertainty 
from one analysis to another, because the analyses han- 
dle jet energies differently, and accept different fractions 
of quark and gluon jets in their respective backgrounds. 
For example, the neural-network jet energy correction 
technique used in the ZH — )- £'^£~bb channels may have 
a different response to the jet energy mismodeling from 
the response in the other H bb channels. Uncertain- 
ties due to initial-state and final-state gluon radiation are 
considered correlated with each other and across chan- 
nels. 



VI. STATISTICAL METHODS 

The results of the searches in each sub-channel are 
represented as distributions of data event counts in in- 
tervals (bins) of a final discriminant variable, which is 
separately optimized for each sub-channel at each value 
of the Higgs boson mass mn. Along with the observed 
data are predictions for each relevant source of back- 
ground, each source of signal, and the associated uncer- 
tainties. These uncertainties affect the predicted yields 
of each component of signal and background, as well as 
the differential distributions of the components in each 
of the histograms. We also consider uncertainties that 
are uncorrelated from one bin to the next of each com- 
ponent of the predictions, usually coming from the lim- 
ited size of Monte Carlo simulated samples. The signal- 
to-background ratios in most of the channels are of the 
order of a few percent or less. The final discriminant 
histograms classify events into categories with different 
signal-to-background ratios. Events with higher discrim- 
inant output values populate bins with larger signal-to- 
background ratios. 

This representation of the search results allows for the 
extraction of constraints on both the signal production 
rates in the decay modes selected and the background 
rates and shapes. Indeed, a major component of the 
sensitivity of the search stems from the ability of the 
data to constrain the rates and shapes of the major back- 
ground sources. The multivariate discriminants sort the 
events based on signal purity. Typically, the low signal- 
to-background portions of the histograms have higher 
statistics and serve to constrain the background rates. 
The shapes of the predictions for each background pro- 
vide the basis by which the background prediction is ex- 
trapolated into the signal-rich region, and the shape un- 
certainties parametrize the extrapolation uncertainties. 

We use the search results to compute upper limits on 
the signal rate for SM Higgs boson production, to deter- 
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mine the best-fit value of the signal strength and cou- 
plings, and to compute p- values for purposes of conduct- 
ing a hypothesis test where the null hypothesis is that a 
Higgs boson signal is absent and the test hypothesis is 
that a SM Higgs boson is present with mass niH- We 
employ both frequentist and Bayesian techniques. The 
upper limits on Higgs boson production and the cross sec- 
tion fits are based on a Bayesian calculation assuming a 
uniform prior probability density of the signal event rate, 



truncated to be non negative. The p- values are computed 
with a frequentist method, although the handling of the 
systematic uncertainties is Bayesian. The approach is 
the same as in Refs. [10, HH- The likelihood function is a 
product over all channels of the Poisson probability of ob- 
serving the data given the predictions, which depend on 
the values of the nuisance parameters that parametrize 
the systematic uncertainties. The likelihood L is shown 
below multiplied by the prior probability density tt, 



L{R, s, b, n, V) X = n n ^^'-1 * ^^^■)"" \ ^ n ^"'"'^^ 
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where the first product is over the number of channels 
{Nc), and the second product is over histogram bins, 
each containing n-ij events. The observed number of 
events in bin j of channel i is n^ . The SM signal pre- 
diction in bin ij is s^, summed over all production and 
decay modes contributing to channel i, and hij is the cor- 
responding background prediction in that bin. The pre- 
dictions Sij and bij are functions of the nuisance param- 
eters V. A nuisance parameter Vk may affect many signal 
and background predictions in a correlated way, such as 
the uncertainty on the luminosity; it may distort the dis- 
tributions of signal and background predictions, as is the 
case with jet energy scales, or it may affect only one bin's 
prediction of one source of signal or background, as is 
the case with Monte Carlo statistical uncertainties. The 
prior probability distributions of the nuisance parameters 
are assumed to be independent Gaussians, and the units 
in which the nuisance parameters are expressed are in 
standard deviations (s.d.) with respect to their nominal 
values. The prior distributions for the nuisance param- 
eters are truncated so that no prediction for the signal 
or background in any channel is negative. The factor R 
is a simultaneous scaling of all signal components. Thus, 
each combination presented here assumes that the rel- 
ative ratios of the contributing Higgs boson production 
and decay modes are as predicted by the model under 
test, within their theoretical uncertainties. We therefore 
present separate combinations assuming the SM, and sev- 
eral choices of models allowing nonstandard couplings. 
We also present separate measurements of R for chan- 
nels that are sensitive to one Higgs boson decay mode at 
a time. 



To calculate the best-fit value of i?, we assume a uni- 
form prior probability density 7r(i?) for positive values of 
R and zero for negative values of R, and integrate the 
likelihood function L multiplied by the prior probabili- 
ties for the nuisance parameters over the values of the 



nuisance parameters 

L'{R,s,b,ri) ^ J L{R,s,b,n,iy)7r{h')dv. (2) 

The best-fit value of R, Rm, is the value that maximizes 
the posterior probability density L'{R)tt{R). The 68% 
credibility interval on R is the shortest interval that con- 
tains 68% of the integral of the posterior density. We 
then define the 95% credibility upper limit on R, Rg^ 
with the following relation: 

0.95= / L'{R)Ti{R)dR. (3) 

We compute the distribution of limits that are expected 
in the hypothesis that no signal is present by simulat- 
ing experimental outcomes and computing i?g5 in each 
of them. The experimental outcomes are simulated by 
varying the values of the nuisance parameters within 
their uncertainties, propagating these to the predictions 
of bij , and then drawing simulated data counts from Pois- 
son distributions with the means of the predicted back- 
grounds. The sensitivity is expressed by the median ex- 
pected limit i?g"°^. A value of Rg^ < 1 indicates that the 
specific signal hypothesis under test is excluded at the 
95% credibility level. 

To evaluate the significance of excess data events com- 
pared with the background prediction, we compute a p- 
value, which is the probability to observe a result that is 
as signal-like or more than the observed result, assum- 
ing that no signal is truly present. A p- value less than 
1.35 X 10""^ is customarily identified as corresponding to a 
three s.d. excess, where the correspondence between the 
p- value and the number of s.d. is computed using the 
integral of one tail of a Gaussian distribution. We rank 
outcomes as more or less signal- or background-like based 
on their i?fit values. We quote the local significance in the 
SM Higgs boson search at mjy = 125 GcV/c^, motivated 
by the recent discovery by ATLAS and CMS 0, S] . 
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FIG. 1: Background-subtracted dijet invariant mass distribution from the combination of all search channels contributing to 
the V Z cross section measurement (left) and collected discriminant histograms, summing bins with similar signal-to-background 
ratio (s/b), for the VZ measurement (right). The expected SM signal contributions are indicated with the filled histograms. 
Normalizations of the subtracted background contributions, with uncertainties indicated by the unfilled histograms, are obtained 
from fits to the data. 



VII. STANDARD MODEL INTERPRETATION 

A. Diboson production 

The search for the Higgs boson at the Tevatron is chal- 
lenging due to large backgrounds relative to the expected 
signal rate. Multivariate techniques are employed to im- 
prove sensitivity and this increases the need to validate 
the background model predictions for rates and kine- 
matic distributions. Over the past few years signals for 
low cross section SM processes have been successfully ex- 
tracted in the same final states as those used for the pri- 
mary Higgs boson searches. For example, the production 
cross section for electroweak single top quark produc- 
tion was measured both in the £vbb 81 1 and J^^bb [82j 
final states which provided important validation for the 
WH -> eM and WH, ZH -> ^^bb searches. Similarly, 
the background model and analysis framework of the 
H W^W~ search have been validated through suc- 
cessful measurement of diboson cross sections have been 
and published in three final states: — >■ W^W~ cross 
section based on the l'^vl~v decay mode [111, jop — » ZZ 
cross section based on the t^l~vv decay mode [8J|, and 
a measurement of the — )- W^Z cross section based on 
the £^i'i~^£~ decay mode [IHl- All three measurements 
were found to be in good agreement with NLO predic- 
tions. 

The searches for WH — > Wbb and ZH Zbb produc- 
tion and decay require careful modeling of large back- 
ground event contributions from Ty-l-jets and Z-|-jets 
production. We gauge the sensitivity of these searches, 
and evaluate the background modeling and analysis tech- 
niques applied within them, by extracting from these 
search channels a combined cross section measurement 
for WZ and ZZ production. The NLO SM cross sec- 



tion for VZ production times the branching fraction for 
Z bb is 0.68 ± 0.05 pb, about six times larger than 
the expected 0.12±0.01 pb cross section times branch- 
ing fraction of VH Vbb for a Higgs boson mass of 
125 GeV/c^. 

This measurement is performed through a combination 
of the same set of search channels used for the WH — > 
evbb, ZH i+£-bb, and WH, ZH ^^bb Higgs boson 
searches. The data sample, event reconstruction, mod- 
cling of signal and background processes, uncertainties, 
and sub-channels are the same as in the Higgs boson 
search. However, dedicated multivariate discriminants 
are trained to separate event contributions of V Z produc- 
tion from those of the other backgrounds and any poten- 
tial contributions from Higgs boson production are not 
considered. Figure [1] shows the background-subtracted, 
reconstructed dijet mass distribution obtained from the 
combination of all search channels. A fit to the data is 
used to determine the absolute normalizations for VZ 
signal and background contributions. 

Separation of the V Z signal component within these 
search channels is obtained from a multivariate discrimi- 
nant that incorporates the dijet invariant mass as one of 
its most powerful kinematic inputs. For improved visu- 
alization of the result of the VZ cross section measure- 
ment, we group event counts from all bins of the final 
discriminant distributions from each of the search chan- 
nels with similar similar signal purity, s/b, and display 
the background-subtracted data contained within each 
grouping as a function of increasing s/b (Fig. [T]). A fit 
to the data is used to determine the absolute normaliza- 
tion of the VZ signal contribution, indicated by the filled 
histogram, as well as the normalization for background 
contributions. The total uncertainty on the background 
prediction is indicated with the unfilled histogram. Based 
on the excess of data events in the highest s/b bins. 
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FIG. 2; Median expected 95% C.L. upper limits on Higgs boson production relative to the SM expectation assuming the 
background-only hypothesis for combinations of search channels within each Higgs boson decay mode and the combination of 
all search channels as a function of Higgs boson mass in the range between 90 and 200 GeV /c? . 



we measure a V Z production cross section of 2.6 'Zit. 
(stat.+syst.) pb, consistent with the SM prediction of 
4.4 ±0.3 pb [zS- 



B. Expected sensitivity to Higgs boson production 

The median expected limit in the absence of signal, 
Rgf^, is shown in Fig. [2] for combinations of the search 
channels within each Higgs boson decay mode, and for 
the full combination of all channels. For Higgs bo- 
son masses below about 130 GeV/c^, searches based 
on the H ^ bb final state provide the greatest sensi- 
tivity. Searches based on _ff — !• W~^W~ arc the most 
sensitive for higher Higgs boson masses. Based on 
the combined result we expect to exclude the regions 
90 < niH < 94 GcV/c2, 96 < mn < 106 GcV/c^, and 
153 < mn < 175 GcV/c^ in the absence of signal. For 
the case of a Higgs boson with a mass of 125 GeV/c^, the 
signal event yields, approximate mass resolutions, and 
median expected limits are shown in Table Hill for combi- 
nations of the channels associated with each Higgs boson 
decay mode. At this mass, H ^ bb has the best sensitiv- 
ity, but the H — )■ W'^W~ searches make an important 



contribution to the combination. 

The final sensitivities of CDF Higgs boson searches are 
a direct result of a substantial effort made over the last 
decade to significantly improve the analysis techniques 
used. The evolution of CDF search sensitivity over time 
is illustrated in Fig. |31 The points show the median ex- 
pected 95% C.L. upper limits on Higgs boson production 
relative to SM expectations assuming the background- 
only hypothesis from the combination of available CDF 
search results performed at various stages over the past 
decade. The integrated luminosities associated with each 
point are the sensitivity- weighted averages of analyzed lu- 
minosities corresponding to the analyzed samples at that 
time. The curves show how the sensitivity of each com- 
bination would be expected to improve in the absence 
of further analysis improvements assuming that sensitiv- 
ity scales inversely with the square root of integrated 
luminosity. With respect to early versions of the CDF 
Higgs boson search sensitivity has been improved by more 
than a factor of two over what would be expected sim- 
ply by incorporating more data. The illustrated gains 
in search sensitivity have originated from a wide array 
of analysis improvements including the inclusion of addi- 
tional triggered events, improved b-]ct identification al- 
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TABLE III: Expected number of signal events, Higgs boson mass resolution, and median expected 95% C.L. upper limits on 
Higgs boson production relative to the SM expectation assuming the background-only hypothesis for combinations of search 
channels within each Higgs boson decay mode at mn = 125 GeV/c^. ^Mass resolution is limited in the H — >■ W'^W~ decay 
mode due to the presence of two neutrinos in the final state, which leads to an under-constrained system. *Mass resolution is 
limited in the ttH — >■ WWbbbb production and decay mode due to the presence of four b quarks in the final state, which leads 
to an ambiguity in jet assignments for reconstructing the Higgs boson mass. 



Channel 


Expected # of signal events 


niH resolution 


Expected limit 
relative to SM 


H ^bb 


87.0 


^15 % 


1.77 




24.2 


Limited^ 


3.25 


_ff — ^ 77 


7.4 


Ri2.5 % 


9.9 


ttH WWbbbb 


3.6 


Limited* 


11.9 


H ^ T+T- 


2.3 


« 25 % 


16.9 


H ^ ZZ 


0.2 


^3 % 


29 




2 4 6 8 10 12 ,14 2 4 6 8 10 12 ,14 

Integrated luminosity (fb ) Integrated luminosity (fb ) 



FIG. 3: Achieved and projected median expected 95% C.L. upper limits on Higgs boson production relative to the SM 
expectation as a function of integrated luminosity, assuming the background-only hypothesis. Each point represents a combi- 
nation of CDF searches performed on the date indicated in the legend. The integrated luminosity associated with each point 
is the sensitivity-weighted average of the analyzed luminosities associated with each contributing channel. The solid lines show 
sensitivity projections, where a scaling inversely proportional to the square root of the integrated luminosity is assumed. The 
information is provided for combined searches performed at mu = 115 GeV/c^ (left) and mu ~ 160 GeV/c^ (right). 



gorithms, implementation of algorithms for improved jet 
energy resolution, inclusion of new search channels such 
as those considering events with additional jets, and im- 
proved multivariate techniques for separating signal and 
background contributions. 



C. Full combination 

The data are categorized into 81 sub-channels for the 
mn = 125 GeV/c^ hypothesis. In order to better vi- 
sualize the results and identify data events causing fluc- 
tuations in the observed limits and p- values with respect 
to expectations for the background-only scenario, wc per- 
form a joint fit of the background predictions for all chan- 



nels to the observed data where nuisance parameters are 
allowed to float within their uncertainties. We then col- 
lect bins from the final discriminant distribution by merg- 
ing bins with similar s/h. The result is shown in Fig. 21 
for the combined channels contributing to the searches fo- 
cusing on the niH = 125 GcV/c^ and 165 GcV/c^ mass 
hypotheses. The predicted Higgs boson contributions 
based on SM expectations summed over the bins with 
similar s/b are shown with the fitted background contri- 
butions overlaid. A subset of the same data is shown in 
Fig. [5] where the data are grouped into wider s/b bins 
and the backgrounds determined from the fit have been 
subtracted. A mild excess of data events is observed in 
the bins with the highest s/h for the niH = 125 GeV/c^ 
hypothesis. No such excess is seen for the case of the 
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FIG. 4: Collected discriminant histograms, summed for bins with similar signal-to-background ratio (s/b), for the combined 
SM Higgs boson searches focusing on the mn ~ 125 GeV/c^ (left) and tuh ~ 165 GeV/c^ (right) hypotheses. Normalizations 
of the background contributions are obtained from fits to the data. Predicted signal contributions, scaled to SM expectations, 
are overlaid on the backgrounds. 
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FIG. 5: Background-subtracted collected discriminant histograms, summed for bins with similar signal-to-background ratio 
(s/b), for the combined SM Higgs boson searches focusing on the mn = 125 GeV/c^ (left) and niH = 165 GeV/c^ (right) 
hypotheses. Background normalizations are obtained from fits to the data, and fit uncertainties are indicated by the unfilled 
histograms. Predicted signal contributions, scaled to SM expectations, are shown with the filled histograms. Uncertainties on 
the data points correspond to the square root of the sum of expected signal and background yields within each bin. 



niH = 165 GeV/c^ hypothesis. 

The likelihood from Equation [T] is used to combine the 
Higgs boson searches from all CDF sub-channels as de- 
scribed in Sec. IVII Figure [5] shows the resulting observed 
upper bound on the signal scale factor Rg^ for potential 
niH values between 90 GeV/c^ and 200 GeV/c^. The me- 
dian expected limit in the presence of no signal, Rgf^, 
is shown by the dark dashed line, while the shaded re- 
gions indicate the limit fluctuation ranges at the level 
of one and two standard deviations. The lighter dashed 



line shows the broad excess in the limits that would be 
expected if a SM Higgs boson with ttih = 125 GeV/c^ 
were present in the data. Values of the observed and 
expected limits are listed in Table IIVI We exclude at 
the 95% credibility level (C.L.) the SM Higgs boson 
within the mass ranges 90 < niH < 102 GeV/c^ and 
149 < tuh < 172 GeV/c^. In the absence of a signal, 
we expect to exclude the regions 90 < m^f < 94 GeV/c^, 
96 < ruH < 106 GeV/c2, and 153 < tuh < 175 GeV/c^. 
Mild excesses in the data compared with fitted back- 
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TABLE IV: Median expected (for the background-only hypothesis) and observed 95% C.L. upper hmits on Higgs boson 
production relative to SM expectations as a function of Higgs boson mass in GcV/c^ for the combination of CDF searches. 



Mass 


90 


95 


100 


105 


110 


115 


120 


125 


130 


135 


140 


145 


Expected 


0.91 


1.01 


0.98 


0.99 


1.06 


1.21 


1.31 


1.46 


1.48 


1.45 


1.35 


1.25 


Observed 


0.45 


0.70 


0.90 


1.12 


1.42 


2.03 


2.82 


2.89 


2.68 


2.22 


2.19 


1.27 


Mass 


150 


155 


160 


165 


170 


175 


180 


185 


190 


195 


200 




Expected 


1.08 


0.94 


0.75 


0.68 


0.82 


0.99 


1.16 


1.49 


1.82 


2.11 


2.37 




Observed 


0.91 


0.71 


0.59 


0.50 


0.85 


1.28 


1.45 


2.31 


3.16 


4.12 


4.79 
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FIG. 6: Observed and expected (median, for the background- 
only hypothesis) 95% C.L. upper limits on SM Higgs boson 
production as a function of the Higgs boson mass for the 
combination of CDF searches. The limits are expressed as 
multiples of the SM prediction for test masses in 5 GeV/c^ 
steps from 90 to 200 GeV/c^ The points are connected with 
straight lines for improved readability. The bands indicate 
the 68% and 95% probability regions where the limits can 
fluctuate, in the absence of signal. The lighter dashed line 
indicates mean expected limits in the presence of a SM Higgs 
boson with mn = 125 GeV/c^. 



ground predictions arc observed, in particular within 
the high s/b bins of the discriminants associated with 
the WH eiybb and ZH i+g-bb searches ^ 
However, in the low mass search region where there is 
overlap with the H ^ bb searches, the H — >■ W~^W^ 
search, which contributes similar search sensitivity at 
TTiH = 125 GcV/c^, docs not contain data excesses in 
the high s/b bins of its discriminants [63l |. By com- 
bining channels, the location of the data excess within 
the range of potential mn values can be partially con- 
strained based on knowledge of the available mass reso- 
lution and expected signal rates from each search chan- 
nel. The constraints are observable in the measured 
values for Rat, which are shown as a function of niH 
along with their associated 68% and 95% C.L. intervals 
in Fig. [71 The moderate excess is localized within the 
region 110 < mn < 140 GeV/c^, where the measured 
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FIG. 7: Best-fit cross section for inclusive Higgs boson pro- 
duction, normalized to the SM expectation, for the combina- 
tion of all CDF search channels as a function of the Higgs bo- 
son mass. The solid line indicates the fitted cross section, and 
the associated shaded regions show the 68% and 95% credibil- 
ity intervals, which include both statistical and systematic un- 
certainties. The mean expected cross section fit values assum- 
ing the presence of a SM Higgs boson at mn = 125 GeV/c^ 
are shown with the dot-dashed line. 



signal rate is found to be consistent with that expected 
from SM Higgs boson production. The best-fit value 
measured for the Higgs boson production cross section 
at TTiH = 125 GeV/c^ is I.54+Q73 (stat.-f syst.) relative 
to the SM prediction. 

The p-valuc is shown as a function of mn in Fig. [51 
The broad excess observed in the cross section measure- 
ment is also visible in the p-value. The p-value for the 
iriH = 125 GeV/c^ hypothesis is 0.0212 corresponding to 
a 2.0 standard deviation excess. A lower p- value (0.0060) 
is observed for the mn — 120 GcV/c^ mass hypothe- 
sis, which is not expected to be distinguishable from the 
rriH = 125 GeV/c^ hypothesis based on the mass resolu- 
tion of the most sensitive search channels. There is also 
approximately a two sigma excess in our data for Higgs 
boson mass hypotheses above «195 GeV/c^. Recent re- 
sults from the LHC experiments strongly exclude the SM 
Higgs boson in this mass range (86l . |87| . Taking this into 
consideration, the mild excess near 200 GeV/c^ is likely 
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FIG. 8: The significance of tlie observed data excess witli re- 
spect to the background-only expectation for the combination 
of all CDF search channels as a function of SM Higgs boson 
mass. The probabilities for the background model to result in 
a best-fit cross section as large or larger than that observed 
in data, p-values, are shown with the solid line. The dashed 
line indicates the mean expected p-values in the presence of 
a SM Higgs boson evaluated separately for each test mass, 
where the associated shaded regions show the ranges of one 
and two standard deviation fluctuations in observed p-values 
for these scenarios. The dot-dashed line indicates mean ex- 
pected p-values for each mass hypothesis in the case of the 
SM Higgs boson with mj/ = 125 GeV/c^. 
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FIG. 9: Summary of best-fit signal cross sections relative 
to SM expectations for the mn ~ 125 GeV/c^ hypothesis. 
Square dots with horizontal uncertainty bars show the fitted 
cross sections obtained from the subsets of CDF search chan- 
nels corresponding to VH — >• Vbb, H — >• W^W~ , H — >■ 77, 
H — >■ T^T^ , and ttH — >■ ttbb production and decay. The solid 
vertical line and associated shaded region illustrate the fitted 
SM cross section obtained from all search channels. 



the result of a statistical fluctuation. 

We study the couplings of a potential SM Higgs bo- 
son by also extracting best-fit signal cross sections for 
different combinations of channels corresponding to spe- 
cific Higgs boson production and decay modes. In par- 
ticular, we perform cross section fits for the subsets 
of CDF search channels corresponding to VH — >• Vbb, 
H W+W-, H ^ J-/, H ^ T+T-, and tiH tthh 
production and decay. Best-fit cross sections relative 
to SM expectations are provided as a function of m^f 
for each of these modes in Table [V] A comparison of 
the individual mode fitted cross sections versus the fit- 
ted SM cross section obtained from all search channels 
is shown in Fig. |9]for the mjy = 125 GcW /c? hypothesis. 
The fitted signal contribution from the H W^W~ 
and H — >■ r+r" channels is zero and for the VH — >■ bb, 
H — > 77, and ttH ttbb channels it exceeds the SM ex- 
pectation. However, all best-fit cross sections are found 
to be consistent within 1.5 standard deviations of SM 
Higgs boson expectations. 



VIII. FERMIOPHOBIC MODEL 
INTERPRETATION 

A number of theoretical models incorporate a Higgs 
boson with couplings to massive bosons as predicted by 
the SM, but negligible or zero couplings to fermions |ll| - 



[l3 |. We denote these as fermiophobic Higgs models 
(FHM). Within these models gg ^ H production is 
negligible, as this mechanism is mediated at lowest or- 
der by quark loops and only higher-order weak interac- 
tions involving W and Z bosons contribute for FHM [s^l . 
Within the FHM interpretation, production rates for 
WH, ZH, and VBF are assumed to be as predicted by 
the SM, while the production rate for ttH is assumed 
to be negligible. Higgs boson decay branching ratios to 
pairs of fermions and pairs of gluons are also set to zero. 
In addition, the decay width r{H — >■ 77) is enhanced 
since quark-loop contributions, which subtract from the 
larger M^-loop contribution, are absent. The complete 
set of decay branching ratios assumed within the FHM 
interpretation arc listed in Table IVll 

Dedicated searches are conducted for — >■ 77 within 
the FHM interpretation to optimize the sensitivity for 
the different event kinematic properties associated with 
the dominant Higgs boson production mechanisms. FHM 
Higgs bosons are produced in association with vector 
bosons, or recoiling from jets in the case of VBF. As a 
result, the Higgs boson px spectrum is shifted to higher 
values for the FHM than the SM, where the dominant 
production mechanism is gg — >■ H. Potential signal con- 
tributions from WH, ZH, and VBF production included 
in the SM i? ^ W+W^ and H ^ ZZ search chan- 
nels are also incorporated. In the H — >■ W~^W~ search 
sub-channel focusing on events with opposite-charge lep- 
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TABLE V: Best-fit signal cross sections, Rat, as a function of tuh for the combination of all SM search channels and for 
combinations of subsets of search channels corresponding to VH — >■ Vbb, H — >■ W^W~ , H — >■ 77, H — >■ t'^t~ , and ttH — >■ ttbb 
production and decay. The quoted uncertainties bound the smallest interval containing 68% of an integral over the posterior 
probability densities, which include both statistical and systematic effects. 
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tons and two or more reconstructed jets, where potential 
signal contributions from these production mechanisms 
are significant, the final discriminant used for the FHM 
interpretation has been re-optimized to focus on the ex- 
pected event kinematic properties of the relevant signal 
processes. The FHM search is performed over the range 
im<mH < 200 GeV/c^ 



No evidence for a fermiophobic Higgs boson is found 
in the data, and upper limits are set on the produc- 
tion rate relative to the FHM expectation. These lim- 
its are shown in Fig. [10] and listed in Table I VII We 
exclude a fermiophobic Higgs boson in the mass range 
100 < ran < 113 GeV/c^, and expect to exclude 100 < 
ran < 122 GeV/c^ in the absence of a Higgs boson signal. 



IX. FOURTH-GENERATION MODEL 
INTERPRETATION AND 
MODEL-INDEPENDENT LIMIT ON gg ^ H 
PRODUCTION 



The lowest-order process mediating the ggH coupling 
in the SM is a quark triangle-loop, with the dominant 
contribution coming from the top quark, and a smaller 
contribution from the bottom quark. The model tested 
here is the standard model with a fourth sequential gen- 
eration of fermions (SM4). The masses of the compo- 
nents of the fourth generation are assumed to be larger 
than the mass bounds from collider experiments. In the 
SM4, the up-type (7/4) and down-type (^4) quarks would 
contribute approximately with the same magnitude as 
the top quark to the ggH coupling, resulting in approxi- 
mately a factor of nine increase in the gg ^ H produc- 
tion cross section and the H ^ gg decay width |15l - [l7| . 
The enhancement is modified by resonant structure in 
the quark loop (the top quark contributes most strongly 
when niH « 2mt), electroweak contributions [l6j . and 
QCD radiative corrections jT^]. The decay branching 
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TABLE VI: Decay branching fractions of the Higgs boson in FHM computed with hdecay Also listed are the observed 
95% credibility level upper limits on the signal rate relative to FHM expectations, and the median expected limits assuming 
no signal is present. 
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fractions of the Higgs boson may further be modified by 
the presence of a fourth neutrino (1^4), which may have 
been too heavy to be discovered at LEP, or due to decays 
to a heavy fourth-generation charged lepton £4. We do 
not include acceptance for these decays in our predicted 
signal yields. The precision elcctroweak constraints that 
place an upper bound on the SM Higgs boson mass [isj 
are significantly relaxed in the SM4, allowing Higgs boson 
masses up to 750 GeV/c^. 

The production cross section for gg ^ H is com- 
puted in Ref. for two scenarios of and rudi, 
but the production rates do not depend significantly on 
these masses, once they are large enough to evade ex- 
perimental bounds. If 2mi^ < mn and 2m,j^ < mn, 
the decay branching ratios have a large impact on our 
ability to test the model. In both scenarios wc assume 
ruu^ = 450 GeV/c2 and ma, = 400 GeV/c^. In the 
first scenario, called the high-mass scenario, wc assume 
rrii^ = mi,^ = 1000 GeV/c^, and in the second sce- 
nario, the low-mass scenario, nii^ = 100 GeV and 
m^^ = 80 GeV/c2. 

We search for gg ^ H production primarily in the 
H — >• W^W~ decay mode, but the H — > ZZ decay mode 
also contributes, particularly for mn > 200 GeV/c^. 
The if — > 77 channels contribute in the SM mainly 
through the gg ^ H production, but this decay mode 



is suppressed due to negative contributions of the quark 
loops relative to the W^-mediated loop in ii 77 de- 
cay. We therefore include only the H W^W~ and 
H — > ZZ — t^£~t^i~ searches in this interpretation. 

The first step is to set a limit on (j{gg — > H) x 
Br{H — VF+VF"), which can be interpreted in a va- 
riety of models. We assume the SM value for the ratio 
of Br{H ZZ)/Br{H W+W-) when combining 
the ZZ results, an assumption which is accurate in the 
SM4. The H — >• W~^W~ channels are re-optimized for 
this search by training the discriminants to separate only 
the gg ~^ H mode from the background, ignoring the 
WH, ZH, and VBF production modes. In setting upper 
limits on the gg ^ H production cross section, we also 
ignore the acceptance for WH, ZH, and VBF produc- 
tion, which yields conservative limits. In setting limits 
on a{gg H) x Br{H W~^W~), we do not include 
uncertainties in theoretical predictions of the production 
cross section or the decay branching ratio, but we include 
the theoretical uncertainties on the relative signal expec- 
tations in the 0-jet, 1-jet, and 2-|~jet event selections in 
the H — > W^W" searches. We search for Higgs bosons 
in the mass range 110 < mn < 300 GeV/c^, in which the 
analysis is expected to be sensitive to the SM4. Limits on 
a{gg ^ H) x Br{H W+W') are listed in Table EIIl 
and arc shown in Fig. 1111 
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TABLE VII: Observed and median expected upper limits on a x Br{H W+W) at the 95% C.L., as well as the predicted 
gg H production cross sections and decay branching fractions in the SM4 with — 80 GeV/c^, me^ — 100 GeV/c^, 
ma^ = 400 GeV/c^ and m^^ = 450 GeV/c^ ' 
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The second step in the SM4 interpretation is to con- 
sider specific model scenarios. In this step we reintro- 
duce the theoretical uncertainties on the predicted cross 
sections due to QCD factorization and renormalization 
scale and PDF uncertainties. The limits obtained are 
shown in Fig. [T^] as multiples of the predictions in the 
two scenarios. In the low-mass scenario, we exclude 
the range 124 < mn < 203 GcV/c^ at the 95% C.L., 
and expect to exclude 123 < mn < 231 GeV/c^. In 
the high-mass scenario, the lack of fourth-generation lep- 
tonic and neutrino decays provides more expected signal 
in the remaining visible decays. We exclude the range 
124 < mn < 206 GeV/c^ at the 95% C.L., and expect 
to exclude the range 123 < mn < 245 GeV/c^. 



X. CONSTRAINTS ON FERMIONIC AND 
BOSONIC COUPLINGS 

Following the recent LHC observations of a new Higgs- 
likc particle with a mass of approximately 125 GcV/c^, 
we focus on this mass hypothesis and test the couplings 
of the new particle, assuming that the mild observed ex- 
cesses in CDF's Higgs boson searches originate from this 
source. Similar studies of the couplings have been per- 
formed by GMS [Hi and ATLAS [8|. 

We assume that the production and decay of the Higgs- 
like particle follows the predictions of the SM Higgs bo- 
son, but with modified coupling strengths to fermions, 
the W boson, and the Z boson. We follow the proce- 
dures and notation of Ref. [oOl , and scale all Higgs bo- 
son couplings to fermions, regardless of flavor, by the 
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FIG. 10: Observed and expected (median, for the 
background-only hypothesis) 95% C.L. upper limits on Higgs 
boson production within the FHM interpretation as a func- 
tion of Higgs boson mass. The limits are expressed as a mul- 
tiple of the expected rate in the FHM for hypothesized test 
masses in 5 GeV/c'^ increments between 100 and 200 GeV/c^. 
The individual points are joined together by straight lines for 
better readability. The shaded bands indicate the 68% and 
95% probability regions in which the hmits are expected to 
fluctuate in the absence of signal. 



factor K/; wc scale the HWW coupling by the factor 
Kyy, and the HZZ coupling by the factor kz- The pre- 
dicted signal rates in each production and decay mode 
are functions of the SM predictions and the factors k/, 
and Kz- The SM predictions are obtained by set- 
ting Kf ^ Kw ^ Kz ~ 1. Because the k factors scale 
the couplings, the production rates and decay widths are 
quadratic functions of the coupling scale factors. The de- 
cay branching ratios are computed from the decay widths 
and thus are ratios of quadratic functions of the coupling 
scale factors. 

For each of the studies described below, we assume a 
uniform prior probability density in one or more of the 
coupling scale factors and compute the posterior proba- 
bility density using all of CDF SM Higgs boson search 
results, integrating over systematic uncertainties. One- 
dimensional intervals arc quoted as the shortest set of 
intervals containing 68% of the integral of the posterior 
density, and the two-dimensional contours are those with 
the smallest areas containing 68% and 95% of the inte- 
gral of the posterior density. The values that maximize 
the posterior probability arc quoted as best-fit values. 

We study both positive and negative values of the cou- 
pling scale factors, although little information on the rela- 
tive signs of the couplings remains after squaring the am- 
plitudes. The posterior probability densities have mul- 
tiple maxima, possibly asymmetric due to interference 
terms in the production and decay in some modes. The 
^77 coupling has a destructive interference term aris- 
ing from the contributions from fermion loops and the 
T4^-boson loop that introduces a term linear in kw and 
Kf. This term breaks the ambiguity of the relative sign 
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FIG. 11: Observed and expected (median, for the 
background-only hypothesis) 95% C.L. upper limits on the 
production rate for gg ^ H ^ W^W^ in picobarns, as 
functions of the Higgs boson mass. The points are joined by 
straight lines for better readability. The bands indicate the 
68% and 95% probability regions where the limits can fluctu- 
ate, in the absence of signal. Also shown are the predictions 
for the two SM4 scenarios described in the text. 
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FIG. 12: Observed and expected (median, for the 
background-only hypothesis) 95% C.L. upper limits on Higgs 
boson production has a function of Higgs boson mass, in the 
SM4 model in the low-mass scenario, which gives the loos- 
est mass bounds. The prediction for the high-mass scenario 
is also shown. The limits are expressed as a multiple of the 
SM4 prediction. The points are joined by straight lines for 
better readability. The bands indicate the 68% and 95% prob- 
ability regions where the limits can fluctuate, in the absence 
of signal. 



between kw and Ky, although the contribution from the 
H — >■ 77 channels is weak in the analyses presented here. 
A smaller interference term exists in the ggH coupling, 
in which the dominant fcrmion-loop contributions inter- 
fere constructively with two-loop clcctroweak contribu- 
tions [13, [11] . A global sign on all Higgs boson couplings 
is unobscrvable in the current analysis. 

We study each coupling scale factor independently, 
holding the others fixed to their SM values, and then 
to study the fits by relaxing the assumptions one at a 
time. We first study kw, setting Kf = kz = \- The 
posterior probability distribution for kw is shown in 
Fig. [131 The factor k\y is constrained to the intervals 
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FIG. 13: Posterior probability distribution for nw from the 
combination of all CDF search channels. In performing this 
fit, the values for nz and k/ are fixed to their SM values 

(kZ = Kf = 1). 



FIG. 14: Posterior probability distribution for nz from the 
combination of all CDF search channels. In performing this 
fit, the values for nw and k/ are fixed to their SM values 

(kW — Kf — 1). 



-1.8 <Kw < -0.8 and 1.0 < kw < 1-7 at the 68% C.L. 
The best fit value for kw is -1.4. We perforin a similar fit 
for Kz, setting k/ = kw = 1. From the posterior prob- 
ability distribution shown in Fig. 1141 kz is constrained 
at the 68% C.L. to the intervals —1.5 < kz < —0.4 and 
0.4 < KZ < 1.5. The best fit value for kz is 1.05. We 
also perform a one-dimensional fit for the Higgs boson 
coupling to fermions, k^, setting kw = kz = 1. The 
posterior probability distribution obtained from the fit 
is shown in Fig. [151 In this case is restricted at 
the 68% C.L. to the intervals -3.8 < k/ < -1.2 and 
2.0 < Kf < 3.0 and has a best fit value of -2.75. 

We also constrain the allowed parameter space in the 
two-dimensional {kw,kz) plane. A fit to the observed 
data is performed allowing all three coupling parame- 
ters to float. Two-dimensional constraints on (kw,kz) 
are obtained from the resulting three-dimensional poste- 
rior probability distribution by integrating over Kf. The 
smallest regions containing 68% and 95% of the integral 
of the posterior probability density are shown in Fig. 1161 
As a result of the global sign ambiguity in the couplings, 
the value of the posterior probability at {—kw,i^z) is 
equal to the value at {kw, —kz). Similarly, the value of 
the posterior probability at {—kw,—kz) is equal to the 
value at {kw, kz). The posterior probability distribution 
in Fig.[Tn]is displayed only for positive values of kw- The 
local maxima within the regions of positive and negative 
Kz are {kw ~ 1.3, kz = 0.9) and [kw = 1-3, kz = -0.9). 

Finally, we constrain the allowed parameter space 
within the two-dimensional (Ky^Kf) plane. Here, kv 
refers to a generic coupling of the Higgs boson to both W 
and Z bosons where the ratio Xwz = / k-z is fixed to 
unity. We compute a two-dimensional posterior probabil- 
ity distribution in the (ty, k/) plane assuming a uniform 
prior probability density. The smallest regions containing 
68% and 95% of the integral of the posterior probability 
density are shown in Fig. 1171 Accounting for the symme- 
tries (— Ky,K/) = (kv,— k/) and {—Kv,—Kf) = {Kv,Kf) 
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FIG. 15: Posterior probability distribution for Kf from the 
combination of all CDF search channels. In performing this 
fit, the values for kw and kz are fixed to their SM values 
{kw — Kz ~ 1). 



the posterior probability distribution is only displayed 
for positive values of kv- The local maxima within the 
regions of positive and negative Kf are {Ky = 1.05, = 
2.6) and {kv = 1.05, Kf = -2.7). 

The results in the (Ky, Kf) plane shown here are more 
constraining than those previously extracted in Ref. [9l| . 
This is due to the inclusion of more channels and the use 
of separate scalings for each signal component, itemized 
by production and decay mode, contributing to individ- 
ual search channels. The search channels with the most 
sensitivity to SM Higgs boson production measure the 
product of vector boson and fcrmion couplings. For ex- 
ample, search modes focusing on decays to fermion pairs 
{bb and r+r^) require production in association with a 
vector boson. Conversely, searches focusing on Higgs bo- 
son decays to W~^W~ and ZZ pairs obtain a majority 
of their sensitivity from gg ^ H production, which pro- 
ceeds mostly via fcrmionic couplings to the Higgs boson. 
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FIG. 16: Two-dimensional constraints in the {kw,k.z) 
plane. The 68% and 95% credibility regions are shown. 
A three-dimensional posterior probability distribution is ob- 
tained from a fit to the observed data in all CDF search 
channels by floating all three coupling parameters {ftw, i^z, 
and Kf). The two-dimensional constraints are obtained by in- 
tegrating the three-dimensional posterior probability density 
over K/. 

Our searches for tiH — > tthh, on the other hand, are sen- 
sitive primarily to k/ in both the production and decay 
modes. Hence, this search channel contributes signifi- 
cantly to the observed constraints on the coupling pa- 
rameters, although it provides only a small contribution 
to combined CDF SM search result. Similarly, search 
channels focusing on events with same-charge di-leptons 
and tri-leptons are sensitive to Ky in both the production 
{VH) and decay {H — )• W^+W^~) modes. These channels 
provide a loose constraint on independent of and 
in the process help eliminate tails in the posterior prob- 
ability distributions. 



XI. SUMMARY 

In summary, we present a final combination of the 
CDF Higgs boson searches. In the context of the stan- 
dard model, we exclude at the 95% C.L. Higgs bosons 
with masses in the ranges 90 < niH < 102 GcV/c^ and 



FIG. 17: Two-dimensional constraints in the (Kv,K,f) 
plane. The 68% and 95% credibility regions are shown. A 
two-dimensional posterior probability distribution is obtained 
from a fit to the observed data in all CDF search channels 
where kv is the generic Higgs boson coupling to W and Z 
bosons assuming the SM value of one for Xwz ~ nw/nz- 

149 < iriH < 172 GcV/c^. In the absence of a signal, 
we expect to exclude the ranges 90 < niH < 94 GeV/c^, 
96 < ruH < 106 GeV/c^ and 153 < mn < 175 GeV/c^. 
An excess of data with respect to background predic- 
tions is observed, which has a local significance of 2.0 
standard deviations at niH = 125 GeV/c^. We ex- 
clude fermiophobic Higgs boson bosons with mass in the 
range 100 < mn < 113 GeV/c^, and expect to exclude 
100 < mn < 122 GeV/c^ in the absence of a Higgs boson 
signal. In the fourth-generation scenario incorporating 
the lowest possible fourth generation lepton and neutrino 
masses, we exclude the range 124 < mn < 203 GeV/c^ 
at the 95% C.L., and expect to exclude 123 < mjy < 
231 GeV/c^. The constraints placed on the fermionic and 
bosonic couplings are consistent with standard model ex- 
pectations. 
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Auxiliary material 
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H->t;V obs 
H-^-cV exp 

ttH^tttb obs 
ttH^tttb exp 

H->Yy obs 
H->YY exp 

VH^Vbb obs 
VH^Vbb exp 
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H^ZZ^4I obs 
H^ZZ^41 exp 

H^W^W" obs 
H^W^W" exp 

Combined obs 
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FIG. 1: Median expected (assuming the background-only hypothesis) and observed 95% C.L. upper limits on Higgs boson 
production relative to the SM expectation for combinations of search channels within each Higgs boson decay mode and the 
combination of all search channels as a function of Higgs boson mass in the range between 90 and 200 GeV/c^. 
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FIG. 2: Background-subtracted dijet invariant mass distribution from the combination of all search channels contributing to 
the VZ cross section measurement. Expected signal contributions from VZ production (red) and a tuh = 125 GeV/c^ SM 
Higgs boson (green) are indicated with the filled histograms. The normalization of the subtracted background contribution, 
with uncertainties indicated by the unfilled histogram, is obtained from a fit to the data. 



